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Abstract. 

We present a general formalism allowing for efficient numerical calculation of the 
production of massless scalar particles from vacuum in a one-dimensional dynamical 
cavity, i.e. the dynamical Casimir effect. By introducing a particular parametrization 
for the time evolution of the field modes inside the cavity we derive a coupled system 
of first-order linear differential equations. The solutions to this system determine 
the number of created particles and can be found by means of numerical methods 
for arbitrary motions of the walls of the cavity. To demonstrate the method which 
accounts for the intermode coupling we investigate the creation of massless scalar 
particles in a one-dimensional vibrating cavity by means of three particular cavity 
motions. We compare the numerical results with analytical predictions as well as a 
different numerical approach. 
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1. Introduction 

The decay of the ground state of quantum field theory, the vacuum, via the creation 
of real particles due to disturbances of the vacuum caused by time dependent external 
conditions demonstrates the highly non-trivial nature of the quantum vacuum. 

The dynamical or non-stationary Casimir effect (see E] and references therein) 
represents one particular example out of the variety of fascinating phenomena occurring 
in the sector of quantum field theory under the infiuence of external conditions [SlIllEl- 
In this scenario the quantum vacuum responds to time-varying boundaries like moving 
mirrors by its decay via the creation of particles out of virtual quantum fiuctuations. 
For calculations involving a single mirror see, e.g., [HI 13 El El • 

In a dynamical (one-dimensional) cavity with one mirror residing at a fixed position, 
say a;i = 0, and a second mirror changing its position X2 according to a given dynamics 
X2it) = l{t) the source of particle creation is twofold. The so-called squeezing of the 
vacuum ilOj due to the dynamical change of the quantization volume (the size of 
the cavity) yields time dependent eigenfrequencies of the field modes. Furthermore, 
boundary conditions imposed on the field inside the cavity at the positions of the 
mirrors cause a time dependent coupling between all field modes. This is denoted 
as the acceleration effect JHl- Being proportional to the velocity l{t) of the boundary 
motion this coupling of the time evolution of all field modes inside a dynamical cavity 
distinguishes the vacuum decay in the dynamical Casimir effect from other quantum 
vacuum effects like particle creation in an expanding universe (see, e.g., [T^) where 
the time dependence of eigenfrequencies is the only source of particle creation. This 
acceleration effect makes it, apart from exceptional cases, impossible to find analytical 
solutions to the field equations. 

A particular example for which analytical solutions are known and the particle 
production can be described by closed form expressions is the scenario of uniformly 
moving boundaries [121 (^^^ ^^^^ El IHI)- Moore |12j found that the mode functions 
(j)k{t,x) inside a one-dimensional dynamical cavity < x < l{t) satisfying both, 
the wave equation [d^ — dl](f)k(t,x) = with 0fc(t,O) = and the non-stationary 
Dirichlet boundary condition 0fc[t,/(t)] = 0, can be written (up to a normalization) 
as 0fe(t,x) = exp{ikTrR(t + x)} — exp{ik7iR{t — x)} provided that R{z) satisfies the 
equation R[t + l(t)] — R[t — l(t)] = 2. This equation can be solved and R{z) evaluated 
exactly, for instance, in the scenario in which l{t) undergoes a uniform motion. 

Another particularly interesting scenario are so-called vibrating cavities |14j . 
Thereby the field is confined between two parallel walls whose distance l{t) changes 
periodically in time: 

l{t) =lo[l + e 6{t)] , with 6{t + T)= 6{t). (1) 

In the case where the frequency of the wall oscillations tUcav is twice the frequency of 
some quantum mode inside the unperturbed cavity, resonant particle (photon) creation 
occurs which makes this scenario the most promising candidate for an experimental 
verification of this pure quantum effect P^ . 
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Let us parametrize the vibrations of the (one-dimensional) cavity as 

Sk{t) = ak sin'' (cut) = sin'' (— (2) 

with Ofc = 1 for k odd and at = 2 for k even. Thereby denotes the frequency of a 
quantum mode inside the unperturbed cavity of length /of- 

For a one-dimensional cavity oscillating with 5{t) = 6i{t) in the main resonance, i.e. 
the frequency of the cavity vibration is twice the frequency of the first field mode inside 
the unperturbed cavity tUcav = 2Q^, analytical solutions for the total particle number 
as well as the number of particles created in the resonance mode k = 1 valid for all 
times have been found in (see also ^H]) under the assumption that the amplitude 
of the oscillations is very small compared to one (e ^ 1). For the same scenario the 
rate of particle creation for higher frequency modes and in the limit ^ 1 has been 
derived in ^Tj by evaluating an approximate solution to Moore's equation. In an 
analytical expression for the number of created particles has been found for more general 
cavity frequencies but short times only, i.e. -C 1. An improved analytical solution to 
Moore's equation is derived in the work [12] where the authors study the energy density 
inside the cavity (see also [201121112211231) and show that the total energy in the cavity 
increases exponentially which was also derived in [T3] . 

The interaction (backreaction) between the cavity motion and the quantum vacuum 
inside the cavity has been studied in, e.g., [211 12S1 120] ■ For work regarding the more 
realistic case of a three-dimensional cavity see, e.g., [101121112111201 ■ 

The analytical results mentioned above have been derived by means of 
approximations such as small amplitudes of the vibrations (e <^ 1) and some of them 
are valid only in particular time ranges. Thus even in the extensively studied vibrating 
cavity scenario there are still open questions. How does the particle creation look when 
the assumption e <^ 1 is no longer valid and the approximative calculations break down 
|. How does particle production behave for higher resonance frequencies as studied 
in but large times? In what manner does the behaviour of the particle production 
change when considering different kinds of cavity vibrations, for instance ^i{t) and 62{t)l 
What happens in the case of detuning (see, e.g., [23(^1011) when the frequency of the 
cavity vibrations does not exactly match the resonance condition? 

In order to answer these questions and to study the vacuum decay in the dynamical 
Casimir effect for a variety of possibly interesting scenarios where less or even nothing 
is known analytically like arbitrary wall motions, cavity vibrations with time-varying 
amplitudes e{t) and massive fields, the problem has to be attacked by using numerical 
methods to account for the coupling of the time dependence of all field modes. 

f This parametrization ensures that the cavity frequency is always twice the frequency of an 
unperturbed mode, i.e. Wcav = 2ri°, and that the total change of the size of the cavity during 
one period is always 2eZo [see also Fig. ^a)]. 

X Note that because of Eqs. and (j^J the amplitude of the velocity of the mirror is ^ oc ewcav and 
thus care has to be taken by increasing the amplitude of the vibrations e such that the velocity of the 
mirror never exceeds the speed of light. 
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Figure 1. Illustration of (a) the cavity vibrations l{t) for fc = 1,2 and 3 studied 
in the present paper and (b) the corresponding function ^[t) — determining the 
time-dependent couphng of all field modes (acceleration effect). 



The goal of the paper at hand is to present a particular parametrization for the time 
evolution of the field modes of a real massless scalar field in a one-dimensional empty 
cavity [0,/(t)] allowing for efficient numerical calculation of the particle production for 
arbitrary motions l{t) of the boundary. A different numerical approach has been recently 
presented in (22] where the author studies the creation of particles in a one- dimensional 
cavity subject to vibrations of the form 62- 

The structure of the paper is as follows. In the next section we briefly review 
the procedure of canonical quantization of a real scalar field in an empty cavity. The 
third section is reserved for the presentation of the formalism yielding a coupled system 
of first-order differential equations from which the number of produced particles can 
be deduced by means of standard numerics. In particular we focus our studies of the 
particle creation in a vibrating cavity on the motions Si(t) and 62(t) (see Fig. Hj). Note 
that the cavity motion Si{t) exhibits a discontinuity in the velocity at the beginning 
of the vibrations, i.e. 6i{t = 0) 7^ and thus l{t = 0) 7^ 0, which can be regarded 
as a pathological feature of the model. However, because of the richness of analytical 
results obtained for this particular scenario (see, e.g., [inillZllIHl) we study the particle 
production for this cavity motion in detail in subsection 4.1 and compare the numerical 
results with the analytical predictions. The impact of the initial discontinuity in the 
velocity of the wall motion is discussed. In 4.2 we study the more realistic cavity motion 
S2{t) with the smooth initial condition Z(t = 0) = 0. We compare the results with the 
scenario discussed in 4.1 and comment on the results and conclusions presented in 
In connection with the results for the cavity motion (52 (t) shown in 4.2 we briefiy discuss 
the effect of detuning in subsection 4.3. Finally we show one example for the cavity 
motion 63(1:) (see Fig. [TJ in 4.4 and conclude in Section 5. 

A more detailed description of the formalism as well as a presentation and discussion 
of numerical results for a wider range of parameters, massive fields and other cavity 
dynamics will be found in [SB]- 
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2. Canonical formulation and quantization 

Let us consider a non-interacting real and massless scalar field ^{t,x) confined to the 
time-dependent interval I{t) = [0,l(t)]. We assume that the scalar field ^(t,x) is 
subject to Dirichlet conditions at the boundaries of I(t), i.e. $(t,0) = = 0. 

The time evolution of ^{t,x) on I{t) is described by the Klein-Gordon equation 
[df — d'^]^{t,x) = which also determines the evolution of the vector potential of 
the electromagnetic field in one space dimension (scalar electrodynamics; see, e.g., 

By introducing an orthonormal and complete set of instantaneous eigenf unctions 
(j)n{t,x) obeying the eigenvalue equation —dl(j)n(t,x) = x) on I(t) we may 

decompose the field as 

= ^g„(t)0„(t,x). (3) 

n 

The eigenfunctions and time- dependent eigenvalues are explicitly given by 

2 



0„(t,x) = W77— 




nir 



nrr 



(4) 



m 

where n = 1,2,... f. Inserting the mode decomposition (jH)) into the Klein-Gordon 
equation yields the set of coupled second-order differential equations [TUl 113 HE] 

+ Mmn{t)qmit) +J2[Mmnit) - Nnm{t)]qm{t) = 0(5) 



m 



where we have defined Nnm = J2k ^nkMmk- The coupling matrices Mnm are determined 
by the integral Mnmif) = /o*"*^ dx4>ni't, x)(f)mit, x) taken over the time-dependent interval 
/(t) (cavity). By using the particular expression for (f)n(t,x) given in ^ one finds 

MUt) = -MUt) = (6) 

l[t) — 

for n m and Mnn{t) = 0. The time evolution of the mode functions depends in two 
different ways on the dynamics of the cavity corresponding to two sources of particle 
creation (TU]: the squeezing of the vacuum due to the non-stationary eigenfrequencies 
(jl)) and the acceleration effect caused by the time-dependent coupling matrix Mnm- 

Quantization is achieved by replacing the classical mode functions by operators, i.e. 
Qn — ^ Qn and Pn ^ Pn, and demanding the usual equal-time commutation relations for 
position and momentum operators. Adopting the Heisenberg-picture the time evolution 
of the operators g„ is described by the same system of coupled differential equations 
(jni) as the classical mode functions. Assuming a static cavity with size Iq for times 
t < the Hamiltonian can be diagonalized by introducing time-independent creation 
and annihilation operators aj^ and a„ of particles with frequency = fi„(t = 0). The 
vacuum state \0,t < 0) = \Qq) which is annihilated by a„, i.e. 

anl^o) = 0, (7) 
is given by the ground state of the diagonalized Hamiltonian H(t < 0) = J2n "-n + ^ 
where the number operator counts the number of particles defined with 

respect to the initial state of the cavity, i.e. for times t < 0. 
f We are using units where Ti ~ c ~ 1. 
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3. Time evolution and particle creation 



We may now use the ansatz 



2Q1 



{t) + a, 



t Am)* 



it) 



(8) 



to parametrize the operator g„ for times t > where the particle operators a„ and 
al defined with respect to the initial vacuum state I^q) have been used as "expansion 
coefficients" [13 ESI- The time dependence is carried by the complex functions e^\t) 
exclusively which obey the same system of coupled second order differential equations 
(0) as the mode functions qn(t). By introducing the functions 



it) 



(m) 



et\t) + Y.Mnkit)et\t) 



it) 



.(m) 



(t) + ^M„,(t)e(r)(t) 



(9) 

(10) 



differentiating them with respect to time t and making use of the second-order 
differential equation (jSJ for e^\t) one obtains the system of coupled linear first-order 
differential equations 

it)= -^[aUt)^t\t)-a-,,itWr\t)_ 

- E + cLit)vt\t)] , (11) 



^(m) 



Vt\t) = - ^ [a~uuit)it\t) - at,it)vt\t) 

-Y.[4nit)^lr\t) + c-,^itWrr\t) 

n 

Thereby we have defined the functions 

'^kit)^^ 



(^kkit) 



1 ± 




and 



Mnkit) ± ^Mknit) 



1 ± 



lit) 



lit) 



\n+k 



n 



k ±n 



(12) 



(13) 



(14) 



for n 7^ /c and c^„(t) = 0. 

Assuming that after a duration ti the motion of the boundary stops and the 
cavity is static again a second set of annihilation and creation operators {An, Al} 
can be introduced to diagonalize the Hamiltonian for t > ti, i.e. H{t > ti) = 
J2n^n ^n^n + \ ■ The grouud state of the Hamiltonian H{t > ti), i.e. the final 
vacuum state = |0,ti), is annihilated by and the number operator iV„ = Al^An 
counts the numbers of particles with frequency fi^ = ^^^(^l) defined with respect to 
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The Bogolubov transformation linking the set of initial state particle operators 
{dn, aj;} with the set of final state particle operators Al} is found to be 



Ar, 



E 



(15) 



where the functions and are linear combinations of ^^\t) and ri^\t) at 

t = ti- In particular 

A^{t,)e:^\t,)+A-{t,)vt\h), (16) 

(17) 



where the function 



A-it) ^ - 



A-(ti)enti)+A+(ii)^f^(ti) 

i{ty 



1 ± 



is somewhat like a measure for the deviation of the final state of the cavity, characterized 
by the cavity length l{ti) = li, from the initial state with Z(0) = Iq. 

Starting from a vacuum state \Qo) the Bogolubov transformation ()15j) has to become 
trivial for times ti = to = 0, i.e., E^^\0) = 2 6nm and H^^'^O) = 0, such that 
^„|f2o) = S„|f2o) = 0, which yields the set of initial conditions 

e)(0) = 25_, vt\0) = (19) 

for the system (fTT|) - (fT^ of coupled differential equations f. 

The number of particles created during the motion of the boundary is given by the 
number of final state particles, counted by Nn = Al^An t, which are contained in the 
initial vacuum state l^o), i-e. 



{no\N^\no) = (u 



(20) 



4 ^ fio 

Knowing the solutions to the coupled system of linear differential equations of first order 
formed by Eqs. (fTT|) and (fT^ this expectation value can be calculated by using the linear 
transformation (fTTj) . 

Accordingly, the total energy E of the created quantum radiation is given by 



no 



(21) 

n=l n=l ^ n=l m=l 

Finding solutions to a coupled system of linear first order differential equations is a 
standard problem in numerical mathematics. By introducing a cut-off quantum number 



and 



t Equations @, 11(J|) and the initial conditions H19|) imply the initial conditions el™''(0) = 
el™'(0) = —ifl^Snrn — Mmn{0) for the complex functions e^\t) which satisfy lO- Therefore, if the 
cavity motion does not start smoothly, i.e. with a non-zero velocity /(O) 7^ 0, such that M„j„(0) 7^ the 
initial conditions for en"'\t) are not simply those of plane waves. This can be seen as well by matching 



quit < 0) = 



/200 



to qn{t > 0) and pn{t < 0) = i 



at e'^'^t 



to 



Pn{t > 0) at i = 0, where q^it > 0) is given by ^ and j5„(t) = q.Jt) + J^m M,i„^{t)q^{t) (see, e.g. 

cm)- 

I Only this particle number operator is physically meaningful for times t>ti. 
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^max the infinite system of coupled differential equations has to be truncated to make it 
suitable for a numerical treatment. Numerical solutions to the remaining finite system 
can be obtained with high accuracy by using standard routines. Finally, the dependence 
of the solutions on the cut-off /Cmax has to be checked in order to guarantee stability of 
the numerical results. In addition, the quality of the numerical solutions can be assessed 
by testing relations like the unitarity of the Bogolubov transformation. Because it is 
beyond the scope of the present paper to go into technical details we refer the reader 
to jnSI for a detailed discussion about the numerics which has been used to obtain the 
results which we show in the next section. 



4. Numerical results 
4.1. 6{t)=6,{t) 

In this section we show some of the numerical results obtained for the scenario of the 
vibrating one-dimensional cavity with l{t) given by 

/(t) = /g [1 + e sin(2fi°t)] (22) 

which has been studied analytically in, for instance, fOl |^ El E] • More numerical 
results are collected in jHS]- We integrate the coupled system formed by the linear 
differential equations (fTT|) and (fT^ for a given interval [0, tmax] and calculate the 
expectation value for all time steps. 

In the authors study the case in which the one- dimensional cavity performs 
oscillations of the form ()22j) with n = 1, i.e. the frequency of the cavity vibrations 
is twice the frequency of the first unperturbed field mode fi^ = ^ inside the cavity. 
Making the assumption of small amplitudes of the oscillations e ^ 1 (trembling cavity) 
they find analytical expressions for the number of created particles in the first mode 
Ni{t) as well as for the total particle number N{t) in terms of complete elliptic integrals 
In particular 

/Vi(() = ^BWA-W-i, (23) 



^2 V 7 V 7 2' 

^ ' 1--k']K'\k)-E{k)K{k) 



where k = ^Jl — exp{— 8r} and t = ^eVl\t = \j^t is the so-called "slow time" 



(24) 



see 



Eq. (6.5) and (6.10) in J3j). These expressions yield N{t) = Ni{t) = for r -C 1 as 
well as N{t) = and Ni{t) = t for r ^ 1. In Figure El (a) we show the numerical 
results for a cavity with initial length Iq = 0.1 and amplitude e = 10~^ obtained for 
an integration time tmax = 3500 and compare them to the analytical expressions ()23j] 
and ()24|) . The numerical results perfectly agree with the analytical expressions of [T3| 
for all times predicting that the initial quadratic increase of both, the total particle 
number and the number of particles created in the resonance mode k = 1, devolves in 
a quadratic increase of the total particle number and a linear behaviour of the number 
of resonance mode particles. The particle spectrum at the end of the integration tmax 
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(a) (b) 




Figure 2. (a) Number of particles produced in a cavity vibrating with (|22|l and n = 1, 
i.e. the main resonance, for parameters Iq — 0.1, e — 0.00001 and /cmax = 15, together 
with the analytical predictions (|23(l and H24|l . (b) Particle spectrum corresponding to 
(a) for the two cut-off parameters fcmax — 10 and fcmax = 15. 



shown in Figure El (b) for the two cut-off parameters fcmax = 10 and 15 indicates the 
stabihty of the numerical results. As one can infer from the spectrum only odd modes 
are created which was also predicted in ^Hj- From the particle spectrum shown for 
short times Nk{t = 250) we read off the value Ni{t = 250) ~ 1.5419 x 10~^ which agrees 
perfectly with the analytical prediction r^|t=25o ~ 1.5421 x 10~^. 

The result of ^H] for r <C 1 has been generalized in ^Hj for more general cavity 
frequencies to 

Nk{T) = {2n- k)kT^ (25) 

for k < 2n and Nk{T) = otherwise, where n characterizes the frequency of the cavity 
vibrations. In Figure El we show the numerical results obtained for the parameters 
of Figure El but cavity frequencies given by = 1.5 and 2 for an integration time 
tmax = 2000. 

The short time particle spectra confirm the analytical predictions. In the first case, 
n = 1.5, the only modes that become excited for times r <^ 1 are the k = 1 and k = 2 
modes and particles of these frequencies are produced in the same amount. From the 
short time spectra shown in Figure El we deduce the numerical values Ni{t = 250) ~ 
3.080 X 10"'^ and N2{t = 250) ~ 3.073 x 10"'^ which are again in perfect agreement with 
the analytical prediction ^ yielding Ni{t = 250) = Niit = 250) ~ 3.084 x lO'^. 

For larger times the behaviour begins to change. The mode k = 1 gets the upper 
hand and slightly more particles in this mode are produced than particles in the mode 
k = 2. Higher frequency modes k > 2 play an inferior role and, at least for this range 
of integration, do not significantly contribute to the total particle number. 

In the second example with n = 2 similar statements hold. As predicted by equation 
()25|) the resonance mode k = n = 2 as well as the modes k = 1 and A; = 3 become excited 
for short times r <^ 1. The amount of particles created in the resonance mode k = 2 
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Figure 3. Number and spectrum of particles created in an oscillating cavity with 
parameters of Figure El but for cavity frequencies given by n = 1.5 and n = 2. The 
particle spectra are shown for two integrations corresponding to two cut-off parameters 
kmax = 20 and 25 for n — 1.5 and fcmax = 25 and 30 for n — 2 to indicate numerical 
stability. 



is slightly larger than the number of particles created in the close-by modes k = 1 and 
k = 3 which are produced in the same amount. As before, this behaviour changes for 
large times. Then the resonance mode k = 2 becomes excited most, followed by the 
close-by modes k = 1 and k = 3, respectively. 

From the numerical simulations (see, e.g., the particle spectra in Fig. ^ we deduce 
that no particles are produced in frequency modes k = 2np with p = 1,2,3, ... where 
n characterizes the frequency of the cavity vibrations. This is a generalization of the 
behaviour found in ^3] that only odd modes become excited in the main resonance 
scenario n = 1 and will be discussed in more detail in 

As next, let us discuss the range of validity of the analytical expressions derived in 
with respect to the assumption e <C 1. For this we consider a cavity with initial 
size Iq = 1 and calculate the number of created particles for amplitudes e covering 
three orders of magnitude. The results of the numerical calculations together with the 
analytical predictions ()23j) and are shown in Figure |3 for e = 0.001, in Figure El 
for e = 0.01 and in Figure (HI for the amplitude e = 0.1§. As one can infer from these 

§ Note that in this rather extreme case with e — 0.1, the maximum velocity of the mirror becomes 
relativistic with - = 2e7r ^ 0.63. 
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Figure 4. (a) Number of particles produced in a cavity vibrating with (|22|l and n = 1, 
i.e. the main resonance, for parameters Iq = I, e ^ 0.001 and fc^ax = 30, together with 
the analytical predictions H23|l and (|24|l . (b) Particle spectrum corresponding to (a) 
for the two cut-off parameters fc,„ax — 20 and fcmax = 30. 




Figure 5. (a) Number of particles produced in a cavity vibrating with (|22|l and n = 1, 
i.e. the main resonance, for parameters Iq = 1, e = 0.01 and fcmax = HO, together with 
the analytical predictions H23|l and (|24|l . (b) Particle spectrum corresponding to (a) 
for the two cut-off parameters fcmax = 100 and fcmax = 110. 



pictures, the rate of particle creation grows rapidly on increasing the amplitude e of the 
oscillations. While for e = 0.001 an integration time of tmax = 500 is needed to obtain 
-^(^max = 500) ~ 0.58, a total particle number of the order of one is already reached 
for if: = 6 in the case of the large amplitude e = 0.1. For large amplitudes like e = 0.01 
and e = 0.1 the number of excited field modes inside the cavity increases drastically. 
This is reflected by the value for the cut-off parameter fcmax which has to be chosen in 
order to obtain numerically stable solutions. Whereas for e = 0.001 the value fcmax = 30 
guarantees stability of the numerical solutions up to tmax = 500 it has to be increased 
to fcmax = 110 to provide stable solutions for e = 0.01 and tmax = 100. In order to 
obtain stable results in the case of the large amplitude e = 0.1 up to tmax = 12 already 
fcmax = 200 modes have to be taken into account for this short integration range. Again, 




Figure 6. (a) Number of particles produced in a cavity vibrating with (|22ll and n ~ 1, 
i.e. the main resonance, for parameters Iq = I, e = 0.1 and /cmax — 200, together with 
the analytical predictions 123|l and (|24|l . (b) Particle spectrum corresponding to (a) 
for the two cut-off parameters fcmax = 190 and fc^ax = 200. 



only odd modes become excited as predicted in [T3] . 

The numerical results for the amplitudes e = 0.001 and e = 0.01 shown in Figure H] 
and Figure respectively, reveal that the expressions ()23|) and derived in by 
means of approximations for e -C 1 describe the numerical solutions very well for all time 
scales under consideration. For e = 0.1 the qualitative behaviour of both, the number of 
particles created in the resonance mode as well as the total particle number, seems still 
to be valid (at least in the shown integration range) but the number of created particles 
is larger than predicted by the analytical expressions. 

For the particular case under consideration it was found that for et ^ 1, i.e. r ^ 1, 
the rate of particle creation in a mode of frequency (k odd) is given by [T7] 

— = -T 26 

dt Tik ^ ^ 

and thus, the number of particles created in the mode k increases linearly for large 
timesf. By expanding Eq. one easily recovers Eq. for the particular case 

k = 1. As mentioned in J3] this asymptotic formula works quite well after r ~ |. 
Because we have already shown that the numerical results for the resonance mode agree 
perfectly with the analytical prediction [see Fig. 121(a) for e = 0.00001, Fig. 111(a) for 
e = 0.001 and El (a) for e = 0.01] we concentrate here on the higher frequencies and 
fig. In Figures m (a) and (b), corresponding to e = 0.01 and 0.1, respectively, we show 
the results for the number of particles created in the modes k = 3 and k = 5 together 
with a linear fit Nk{t) = akt + bk to the numerical values for certain time ranges. For 
the amplitude e = 0.01, for which the slow time value r = 1 corresponds to t ~ 62, the 
rate of particle creation obtained by fitting the data for times r > 1 agrees very well 
with the values predicted by Eq. ()26|) as one infers from Fig. [71 (a). From our numerical 
calculations we find the values 03 = 0.00417 and 05 = 0.00254 which are in very good 

f Note that in XI-. ^ factor 2 was missed which has been corrected in |15| . 
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Figure 7. Number of particles created in the modes fc = 3 and fc = 5 for r > 1 and 
parameters = 1 and (a) e = 0.01 and (b) e = 0.1, corresponding to Fig. [Hl(a) and Fig. 
ini(a), respectively, together with a fit of the data to the linear law Nk{t) — att + bk- 



agreement with the values = ^ ^ 0.00424 and 05 = ~ 0.00255 predicted by Eq. 
(I26|) . The integration range used in the numerical simulation for e = 0.1 yielding the 
results shown in Fig. El allows us as well to compare the data to the prediction (|26|) 
for times r > 1. In this case, where the assumption e -C 1 is no longer valid, we still 
find a relatively good agreement of the numerically obtained rate of particle creation 
with the predicted one [ see Fig. [7| (b) \]. Interestingly, whereas the numerical result 
for the number of particles created in the resonance mode is not very well described 
by the analytical prediction ()23p for times r > 1 [see Fig. IHl (a) which shows that 
the rate of particle creation is larger than the (in the limit e ^ 1) predicted one] the 
numerical results for the higher frequency modes k = 3 and k = 5 agree comparatively 
well with Eq. The range of integration used in the simulations for e = 0.00001 

and e = 0.001 (see Figs. ElandlH) is not large enough to compare the numerical results 
with the analytical predictions for times r > 1. For the amplitude e = 0.00001 [see Fig. 
121 (a)] the number of particles created in the modes k = 3 and = 5 is still in a phase 
of acceleration up to the maximum integration time tmax = 3500 which corresponds to 
r ~ 0.55. Similar statements hold for the numerical results obtained for e = 0.001 [see 
Fig. m (a)] where the maximum integration time corresponds to r ~ 0.79. However, as 
discussed before, we have found a perfect agreement of the numerical results obtained 
for e = 0.01 with the analytical expression for times r > 1 [see Fig. (3 (a)]. Even 
for the large amplitude e = 0.1 the numerical results coincide with the predictions quite 
well [see Fig. [71 (b)]. Therefore, our numerical simulations show that the rate of particle 
creation for higher frequency modes is very well described by the analytical expression 
(Uni) beginning at r ~ 1 f. 

Now, let us have a more detailed look at the process of particle creation. In Figures 
EJa) - Efa) we have included additional pictures showing the total particle number and 

I The scattering of the numerical values in this case is due to the oscillations in the particle number 
which can be also seen in Fig. [HI (a) . 

f It was noted in ^21 that Eq. H26(l is valid only for not very large numbers k due to limitations of the 
used approximations. For e — 0.01 we have found that Eq. perfectly describes the rate of particle 
creation for fc = 7 and fc = 9 as well. 
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the number of created resonance mode particles for short times in high time resolution. 
For e = 0.01 and 0.1 we have illustrated the background dynamics l{t) as well. For 
e = 0.001 [see Fig. El^a)], again, the agreement of the numerical results with the 
analytical prediction N{t) = is convincingly illustrated. In these high resolution 
pictures oscillations in the particle number correlated with the motion of the boundary 
become visible for e = 0.01 and 0.1 J. This observation relies on the fact that in the 
numerical calculation we evaluate the expectation value ()20|) at every time step of the 
integration and not only for times at which the instantaneous state of the cavity equals 
the initial one §. Clearly, the larger the amplitude of the cavity vibrations, the larger 
the amplitude of the oscillations in the particle number. For e = 0.01 the numerical 
result for the number of particles created in the resonance mode agrees well with the 
analytical expression (j23|l but the numerical values for the total particle number exceed 
the analytical prediction. Similar statements hold in the case of the large amplitude 
e = 0.1 [see Fig. Efa)]. In both cases one observes a jump in the total particle number 
from zero to a much larger value at the first time step of the integration. Consequently, 
field modes of higher frequencies k = 3,5,.. become excited even at the first step 
of integration and contribute to the total particle number from the very beginning 
(remember that for e = 0.1, k^ax = 200 modes had to be taken into account to provide 
numerical stability for tma.^ = 12). 

The excitation of high frequency modes from the very beginning may be due to 
the fact that the cavity motion (j^^ is not smooth at t = but starts with a non-zero 
velocity which is proportional to the amplitude e ||. The discontinuity in the velocity 
of the mirror at t = then induces the excitation of modes of higher frequencies and 
therefore acts as a source for spurious particle creation which manifests itself in a kind 
of particle background which is present right from the beginning. This effect is nicely 
illustrated in Figures|Hl(a) and|ni(a) which show that the total particle number oscillates 
on top of this background. Being proportional to e, this effect should play a secondary 
role for tiny amplitudes e -C 1 but become more and more important as e increases. 
This is confirmed by our numerical simulations showing that for tiny amplitudes like 
e = 0.00001, the total particle number for short times is practically identical to the 
number of particles created in the resonance mode and therefore the excitation of higher 
frequency modes does not appear. In contrast, instantaneous particle creation takes 
place at the first integration step for the amplitudes e = 0.01 and e = 0.1 [see Figs. El 
(a) and|Bl(a)] due to contributions from modes of higher frequencies to the total particle 

I For e ~ 0.001 oscillations in N{t) appear as well but with a rather tiny amplitude, not visible in 
Figure ^a). 

§ If we ask for the number of created particles at times t = nT only, where n is the number of cavity 
oscillations of period T, we can use the function 77^™'' directly to evaluate the expectation value H20|l . 
because A+(nT) = 1 and A-{nT) = [see Eqs. ^ and lfTH|l ]. 

II In ^112] it was found, that due to such an initial discontinuity in the velocity of the wall the energy 
density inside a vibrating cavity develops d function singularities. Furthermore, the discontinuous 
change in the velocity of the uniformly moving mirror discussed in |12| leads to a logarithmically 
divergent particle number (see also |35p. 
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number, induced by the discontinuity in the velocity of the mirror at t = 0. 

Let us summarize the resuhs obtained for the cavity motion (j221)- For this purpose, 
the numerical results for the main resonance scenario, i.e. cUcav = ^Q^, have been 
arranged in Figure |H1 (a) together with the analytical predictions (j^Hj) and ()24|) . The 
numerical calculations reveal that the analytical expressions ()2H|1 and fl^i describe the 
numerical results perfectly for the amplitudes e = 0.00001 and e = 0.001 for all time 
scales under consideration. Furthermore, for e = 0.01 where in high time resolution the 
effect of the instantaneous particle creation at the first integration step becomes visible 
[see the high resolution picture in Figure|Hl(a)] the behaviour of the total particle number 
as well as the number of particles created in the resonance mode is still in very good 
agreement with the analytical predictions. For the certainly somewhat artificial scenario 
with e = 0.1 where the mirror starts its oscillations instantaneously with a relativistic 
velocity (recall that in this case w/c ~ 0.63) the analytical expressions and do 
not match the numerical results. 

In [in] the authors found, in addition to the analytical expressions and (j^ . 
the closed form expression 



for the energy of the created quantum radiation which grows much faster than the 
total particle number. For a calculation of the energy density inside a vibrating cavity 
see, e.g., [121 HH 123 • In Figure |Sl(b) we compare the analytical prediction with 
the energy of the created quantum radiation calculated numerically by means of Eq. 
()21|) . For the small amplitudes e = 0.00001 and e = 0.001 the numerical results agree 
perfectly with the analytical prediction (jTTjl . In the case of e = 0.01 the numerically 
calculated energy of the created quantum radiation deviates slightly from the analytical 
prediction for small as well as large times where the total particle number is still in 
very good agreement with the analytical expression [cf. Figs. Efa) and IHfa)] 
The energy exceeds the predicted value for small times indicating again that in this 
case the effect of instantaneous excitation of higher frequency modes due to the non- 
smooth beginning of the cavity motion (j^^ seems to become important. For e = 0.1 the 
numerically calculated energy of the created quantum radiation deviates drastically from 
the analytical prediction for short times showing the effect of instantaneous excitation of 
high frequency modes caused by the discontinuity in the velocity of the mirror at t = 
quite impressively. Right from the very beginning much of the energy of the cavity 
motion is transfered to quantum modes of higher frequencies yielding a comparatively 
large energy of the quantum radiation even for small times. 

Let us now proceed to the investigation of a (more realistic) scenario in which the 
vibrations of the cavity start with zero velocity. 

^ Note that the energy of the created quantum radiation (|21|l is much more sensitive to contributions 
from higher frequency modes than the total particle number N = J2n ^n, due to the multiplication of 
the number of particles created in the mode n with the frequency . 




(27) 
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Figure 8. (a) Summary of numerical results obtained for the particle production in 
a cavity oscillating with H22|l and n = 1, i.e. the main resonance. We show the total 
particle number as well as the number of particles created in the resonance mode for 
amplitudes e = 0.00001 for an initial cavity size Iq = 0.1 corresponding to Fig. |21 
and e — 0.001, 0.01 and 0.1 for Iq = 1, corresponding to Figs. El and compare 
them with the analytical expressions H23|l and (|24(l . (b) Comparison of the numerical 
results obtained for the energy associated with the created quantum radiation for the 
parameters shown in part (a), with the analytical expression (|27|l . 



12. 6{t)=62{t) 

In this section we study the cavity motion 



l(t) =lo\l + 2e sm^{nl t)] = /q [l + e (1 - cos(2 t) 



(28) 



with smooth initial conditions Z(0) = Iq and /(O) = which has also been studied 
numerically in In the following we compare the process of particle creation in 

a one-dimensional cavity vibrating with ()28j) with the results obtained for the cavity 
motion (|^. Therefore we set /q = 1 and restrict ourselves to amplitudes e = 0.001, 
0.01 and 0.1. The results of the numerical simulations are shown in FigureslHl-ITT] where, 
for reasons of comparison, we have included the analytical expression (j^^ as well. 

We observe that the process particle creation caused by the cavity motion (j28| is 
different from the one driven by (j22I)- The qualitative behaviour of the total particle 
number shows no differences for relatively short times compared to the background 
motion (j22I), i-e. N{t) grows quadratically with time. But for larger times this behaviour 
begins to change. In the case of e = 0.001 [see Figure|ni(a)] the initial quadratic increase 
of the particle number is followed by a slowing down of the rate of particle creation with 
time and the particle number enters a region in which it is effectively described by a 
linear behaviour. Thereby the total particle number is always less than the number 
of resonance mode particles created in a cavity vibrating with ((221) • For e = 0.01 [see 
Figure Uni (a)] and e = 0.1 [see Figure ITT] fa)] the same qualitative behaviour is observed. 
But for larger times |, the particle number seems to have a tendency to leave the linear 

I Note that the expressions "short times" and "large times" extensively used throughout the paper 
refer to the so called " slow time" r = ^ e 0° i as introduced in P31 rather than to the time variable t 
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Figure 9. (a) Number of particles produced in a cavity vibrating with H28|l and n — 1 
for parameters Iq — I, e — 0.001 and fcniax — 30, together with the analytical prediction 
l)23|l . (b) Particle spectrum corresponding to (a) for the two cut-off parameters 

^max = 20 and fcmax = 30. 
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Figure 10. (a) Number of particles produced in a cavity vibrating with (|28ll and n = 1 
for parameters Iq = 1, e = 0.01 and /cmax = 50, together with the analytical prediction 
(|^ . (b) Particle spectrum corresponding to (a) for the two cut-off parameters 
kmaK = 40 and fc^ax = 50. 



regime due to a further deceleration in the rate of particle production. The number of 
particles created in the mode k = 1 shows this slowing down very clearly. 

We have found that the qualitative behaviour of the total particle number can be 
described very well by a linear law in the time range 1/2 < r < 1. This is illustrated in 
Fig. nowhere we show fits of the numerical results to N{t) = a t + b in the corresponding 
time ranges. From these plots the linear behaviour of the total particle number in the 
time range ^ < t < 1 becomes evident §. 

used in the numerical simulations. Roughly spoken, we use the term "short time" ("large time") for 
T < 1 (r > 1). 

§ From the parameters a obtained by fitting the numerical results in the time range 1/2 < r < 1 one 
may deduce, at least as a good approximation, the dependence a ^ e. Later on we will see that, more 




Figure 11. (a) Number of particles produced in a cavity vibrating with (|28|l and n ~ 1 
for parameters = 1 ,e = 0.1 and fc^ax = 80, together with the analytical prediction 
H23|l . (b) Particle spectrum corresponding to (a) for the two cut-off parameters 

A:max = 70 and fcmax = 80. 
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Figure 12. Total particle number in the time range ^ < t < 1 for parameters Iq = 1 
and (a) e = 0.001, (b) e = 0.01 and (c) e = 0.1, corresponding to Figs. IKa) -ITTlfa). 
together with a fit of the data to the linear law N{t) — at + b. 



To demonstrate the transition from the quadratic to the hnear behaviour of the 
total particle number more carefully we show a summary of the numerical results for 
the total particle number in Fig. (a) together with the analytical prediction ()24p 
[derived for the cavity motion ((221)] as well as linear functions N{t) = at + b with 
parameters a and b corresponding to Fig. El As on infers from Fig. El (a), the total 
particle number is very well described by Eq. up to t ~ 100 for e = 0.001 and up to 
t ~ 10 for e = 0.01. Then the deceleration in the rate of particle creation sets in and the 
time evolution of the total particle number enters the linear regime 1/2 < r < 1. While 
for e = 0.001 the particle number stays in the linear regime by reaching the end of the 
shown integration range the further deceleration of the rate of particle creation becomes 
visible for e = 0.01 and 0.1, and the particle number deviates more and more from the 
linear behaviour by approaching the maximum integration time. To investigate this 
deceleration in the rate of particle creation in more detail, one clearly has to consider 

generally, cl — seems to hold such that, in the linear regime, ^^^j^-^ = § (for n = 1, i.e. Wcav = 2D,i). 
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larger integration times. We will address this question later on. 

The particle spectra shown in Figures IHl (b) - ^2 (b) reveal that the characteristic 
features remain unchanged compared to the cavity motion (j22I), e.g., only odd modes 
are created. But the number of particles produced in each mode is smaller compared to 
Figures m (b) -ini(b) yielding a smaller total particle number as described above. While 
for e = 0.001 the same cut off parameter k^^x has to be used for both cavity dynamics, 
the number of modes which have to be taken into account to ensure numerical stability 
for amplitudes e = 0.01 and 0.1 is smaller for the cavity motion ()28|1 compared to (j221), 
indicating that fewer modes of higher frequencies are excited in a cavity oscillating with 
(PH|l compared to cavity vibrations of the form \. 

In Figures El (a) - ^2 (a) we have included pictures showing the numerical results 
for short times and with a high time resolution which we now compare with the 
corresponding results obtained for the cavity motion ()22j) shown in Figures E] (a) - El 
(a). By comparing the high resolution picture in Fig. |H1 (a) with the one in Fig. HI 
(a) we clearly see (in this resolution) no difference. For both cavity motions, the total 
number of created particles is well described for short times by N{t) = r^. For e = 0.01 
and e = 0.1 the oscillations in the particle number, caused by the fact that we evaluate 
the particle number at every integration step as explained in the former subsection, now 
mimic the background motion ()28|). They start smoothly and differ from the oscillations 
in a cavity driven by (j^^ only by a phase. Apart from that, the number of particles 
created in the mode of frequency behaves exactly the same in both cases. For the 
total particle number the situation is now different. Without the discontinuity in the 
velocity of the mirror at the beginning of the integration the total particle number does 
not exhibit the jump at the first step of integration, which is characteristic in the case 
of the background motion (|22I)- By comparing the short time pictures in Fig. El (a) 
and Fig. ^1 (a) for e = 0.01 as well as in Fig. El (a) and Fig. ^2 (a) ioi e = 0.1 we 
clearly recognize the effect of spurious particle creation caused by the discontinuity in 
the velocity of the mirror motion ()22j) at t = 0. The amplitudes of the oscillations in 
the total particle number itself are of the same height but for the cavity motion ()22|1 the 
oscillations sit on top of a kind of particle background which is present from the first 
step of integration as described in the last subsection. For the background motion p8|) 
with the smooth initial condition /(O) = no instantaneous particle creation takes place, 
the particle background does not appear and thus, the total particle number starts to 
increase smoothly during the first integration steps. 

In Figure (131 (b) we show the energy associated with the produced quantum 
radiation corresponding to Figures El - (HI calculated by means of Eq. (j2T| together 
with fits of the numerical data to a power law E{t) = (3t". We find that in the time 
ranges under consideration the energy E{t) produced in a cavity subject to vibrations 

I Note that the cavity motion ^ [and therefore H28I) ] was chosen such that the total change Al of 
the cavity length is 2e for all k to have comparable situations. See also Figure^ If Al were different, 
Al — 2e for H22I) and Al = e for (|28|l . say, the total number of particles created in both cases would 
differ drastically (up to an order of magnitude) as we have observed in our simulations 




Figure 13. (a) Summary of the nmnerical results obtained for tlie particle production 
in a cavity oscillating with (|28|l and n — 1. We show the total particle number for 
Iq = 1 and amphtudes e = 0.001, 0.01 and 0.1, corresponding to Figures. EKa) -ITTIfal 
and compare them for short times with the analytical expression H24|l and for larger 
times with the linear law N{t) = at + b with parameters a and b given in Figures IT^fal 
- (c). (b) Time evolution of the energy associated with the created quantum radiation 
corresponding to the parameters shown in part (a) together with a fit of the numerical 
data to the power law E{t) = fif^. 



of the form ()28|1 and amplitudes e = 0.001 and 0.01 is very well fitted by a = 2 
f. Thus the energy increases quadratically for e = 0.001 and 0.01 in contrast to the 
exponential behaviour of the energy for the cavity motion [see also Eq. (P7j) and 
Fig. IHl (b)]. For the amplitude e = 0.1 the rate of energy production is even less (in the 
discussed time range) and described by the exponent a = 1.7346 \. Therefore, the rate 
of energy production is much smaller in a cavity vibrating with ()28j] compared to (j221)- 
Partly, this is due to the fact that, as mentioned above, fewer modes are excited in the 
case of the background motion (j^HJ. Furthermore, as one can infer from the presented 
particle spectra Nkit^ax), the number of particles A^^ created in each excited mode with 
frequency at a given time in a cavity vibrating with ()28|) is less compared to the 
corresponding number of particles created in a cavity subject to the background motion 
(El §• 

Now let us compare our method with a completely different numerical approach. 
In the author studies the creation of massless scalar particles in a one- dimensional 
cavity subject to the motion (pSj) by solving the Klein - Gordon equation using an 
improved Leap Frog algorithm (consult the paper for details). Figure IT^ summarizes 
our numerical results obtained for n = 1, i.e. c^cav = 2^, and the parameters Iq = 50 
and e = 0.02 which correspond to the parameters used in the simulations of [S2]- We 

f For e ~ 0.001 and 0.01 the parameters a and /3 were obtained by fitting the power law to the 
numerical data in the range imax/10 < t < fmax- 

I In this case the fit to the numerical data was done in the range 5 < < < imax = 12. 

§ This statement does not hold, of course, for short times for which the particle number behaves 

qualitatively similar for both cavity motions. 




Figure 14. (a) Comparison of tlie particle production for tlie two cavity motions H22(l 
and (|28|l with n — I, i.e. the main resonance, and parameters Iq = 50 and e = 0.02 
corresponding to the parameters studied in [IV2\ . The analytical predictions (|23|l and 
(|24|l . valid for the background motion H22|l . are shown as well, (b) Particle spectra 
corresponding to (a). 



performed our calculation for an integration time tmax = 1400 which is twice the one 
considered in [32] ■ For a better comparison, we depict our numerical results for the 
shorter integration range [0 < t < 700] in an additional picture in Fig. El (a) which 
should be compared with Fig. 5 in [221 and demonstrates that our results perfectly 
agree with the results presented in |^. Furthermore, the outcome of the numerical 
simulations obtained for the same parameters but for the cavity motion (j^^ is shown 
in Figure [T3I (a) as well, which perfectly agrees with the analytical expressions and 
(1^ also indicated in this plot. 

The result for total number of created particles shown in Fig. 5 of for an 
integration time tmax = 700 is fitted by a power law oc with a = 1.85 (in the range 
100 < t < 700) which is interpreted by the author as in reasonable agreement with 
the prediction a = 2 derived in [T3] for the cavity dynamics (|^. but not (PHj) ! As 
we have shown with our numerical simulations, the two cavity motions (|^^ and (j2Hl), 
apart from comparatively short times, do not yield the same behaviour for the number 
of created particles §. As one infers from Figure El (a) the behaviour of the total 
number of created particles is very well described by Eq. (j^ . i.e. a quadratic increase, 
up to t ~ 300. After that, the rate of particle creation slows down and the particle 
number shows the transition to the linear regime as discussed before in detail. Thus, 
the exponent a = 1.85 found in |32j by fitting the total particle number to a power law 
for 100 < t < 700 is explained by the fact, that this time range contains the transition 
from the initial quadratic increase of the particle number to a linear behaviour which 
becomes visible only for larger integration times as confirmed by our simulations [see 
Figure IT¥a)]. 

§ Recall that the analytical expression H24|l for the total number of created particles derived for the 
cavity motion (|22|l behaves like t'^ for small as well as large times. 




Figure 15. (a) Total number of particles produced in a cavity vibrating with l|28(l 
and n = 3, i.e. oJc&v = 2 fig, for parameters = 50, e = 0.02 and /cn^ax = 80. (b) 
Particle spectrum corresponding to (a) for the two cut-off parameters fcmax = 70 and 



Nevertheless, the work [22] provides a good possibihty for checking our method 
against a completely different numerical approach. To illustrate this by means of 
one more example we show the numerical result for = 3 , i.e. cUcav = 2 fig (which 
corresponds to the case = 6 in [S21)! and the same parameters Iq = 50 and e = 0.02 in 
Figure El which perfectly agrees with the corresponding result obtained in As 
mentioned by the author of this behaviour of the particle number does not fit any 
simple expression. 

Now, let us finally address the issue of further deceleration of the rate of particle 
creation which yields a deviation of the time evolution of the number of created particles 
from the linear behaviour for larger times which is indicated in Figs. ^1 (a) and ^2 (a) 
[see also Fig. El (a)]. For this purpose it turns out that, from a numerical point of view, 
it is convenient to maintain the parameters Iq = 50 and e = 0.02. For these parameters. 
Fig (a) nicely demonstrates the transition from the initially quadratic to the linear 
behaviour of the number of particles created in a cavity vibrating with (j^H|) and n = 1, 
i.e. uicstv = In Figure IT^ we show the numerical results obtained for an integration 
time tmax = 4000 (which corresponds to Tmax = > 1) together with Eq. and a 
linear fit to the total particle number f. As already observed in Fig. El (a) the quadratic 
increase of the total particle number up to t ~ 300 is followed by a transition to the 
linear regime. For times t > 1600 (r > 1) the particle number starts to leave the region 
in which it is effectively described by a linear behaviour and shortly after (t ~ 2000) 

t Compare Figure [Tsl (a) with the case n = 6 in Fig. 6 of [32 and Figure ^(b) with Fig. 7 of jSH- 
f The time interval in which we have performed the fitting procedure is 796 < t < 1592 corresponding 
to the "definition" of the linear regime 1/2 < r < 1 given above. The value a = 0.000385 0.0004 
may indicate the validity of the dependence a = j-, i.e. "^^^^"^ = ^- We have performed numerical 
simulations for Iq = 50, e = 0.03 and = 25, e o!o2 up to r = 1 yielding a = 0.00058 - 0.0006 and 
a = 0.00077 ^ 0.0008, respectively. These results strongly support the indication that = ^ holds 

in general for 1/2 < r < 1. 
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Figure 16. (a) Number of particles produced in a cavity vibrating with (|28|) and 
n = 1 for the parameters of Fig. but now for an integration time imax = 4000 
and fcmax = 70. The analytical expression H24(l and the linear law N{t) — at + b with 
a = 0.000385 obtained by fitting the numerical data in the time range | < r < 1 are 
shown as well (see the discussion in the text), (b) Particle spectrum corresponding to 
(a) for the two cut-off parameters fcmax = 60 and fcmax = 70. 



the further deceleration of the rate of particle creation yields a rising deviation of the 
numerical values from the linear fit. The number of particles produced in the modes 
k = 1 and k = 3, shown in Fig. (a) as well, exhibit the same qualitative behaviour. 
Therefore, for a cavity vibrating with ()28j) . the overall time evolution of the total particle 
number as well as the number particles created in the mode k = 1 cannot be described 
(and fitted) by a comparably simple expression similar to the analytical results (f^H|) and 
fl24|) valid for the cavity motion ()22|) . This behaviour is due to a detuning effect as we 
will explain in the next subsection. 

4-3. 6{t) = S2{t) - Detuning and the "true resonance" 

In order to understand the long time behaviour of the number of particles produced 
in a cavity oscillating with ()28|1 we have to formulate the resonance condition more 
carefully. For a resonance to happen the external frequency Wcav has to be twice the 
eigenfrequency of a quantum mode defined with respect to the average position / of the 
cavity, i.e. c^cav = 20„ where Cln = nir/l. 

Whereas for the cavity motion ()22|) the average frequency Qn corresponds to the 
frequency n° defined with respect to the initial vacuum state \Qo) (see section 2) the 
situation is different for the background motion (j2HI)- In this case the average position 
of the cavity is / = /o(l + e) (see also Fig. [T)) and therefore n„ 7^ Thus, Eq. (j2HI) 
is not a real resonant cavity motion because the cavity frequency does not match the 
exact resonance condition. Such an effect is called detuning (see, e.g., [23110]) and can 
be parametrized by the detuning parameter A. In particular, for the cavity motion ()28j] 
we have cUcav = 2(n„ + A) with A = fl^ — fin = f^°(^— lo)/l, which for = 1 and Iq = I 
reduces to A = tt ( -rr- ) ~ vr e. 
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We have studied how detuning affects the particle production in a one- dimensional 
cavity in dependence on the detuning parameter A in detail. For instance for the cavity 
motion (j22I), we have found that the number of created particles oscillates in time 
with a period and an amplitude depending on the detuning parameter A, i.e. particles 
are created and annihilated periodically in time. Therefore, detuning causes phases of 
acceleration as well as deceleration in the time evolution of the particle number and is 
the reason for the qualitative behaviour of the particle number which we have observed 
in the former subsection. These results will be reported and discussed elsewhere in more 
detail 

In this paper we restrict ourselves to the demonstration, that on replacing the cavity 
frequency tUcav = = 2^ in (PHjl by c^cav = 2Qi = 2 i-e. A = (no detuning), 

the two cavity motions (j22|) and (j28|) do indeed yield the same qualitative long time 
behaviour for the particle number f. In Figs. El and ^1 we show the numerical results 
for the parameters /q = 1 and e = 0.001 and 0.01, respectively. For a better comparison, 
we have depicted the analytical expressions and derived for the cavity motion 
(j^^ as well. Thereby we now use r = ^eClit as definition for the "slow time" to account 
for the modified cavity frequency. 

As we infer from Figs. El (a) and El (a) the time evolution of the number of 
created particles is now qualitatively as well as quantitatively in good agreement with 
the analytical predictions (j23|) and (j24j) . Instead of the linear behaviour and the further 
deceleration in the rate of particle creation observed for the cavity motion (j^Hj) with 
i^cav = 2fl1 [cf. Figs. M (a) and |TU| (a)] the time evolution of the number of created 
particles shows now the typical characteristics of resonant particle creation. For this 
reason we call the scenario with the cavity dynamics (j28|l and Ucav = 2fli the "true 
resonance" scenario. 

Note that the number of excited modes (the value of the cut-off parameter fcmax) 
is now equal to the number of modes excited in a cavity vibrating with ()22j) [compare 
Figs. HHI (b) and El (b)]. Thus, the fact that the number of excited modes is less for the 
cavity motion with cUcav = 20? is due to the effect of detuning and not because of 
the smooth initial condition. This indicates that the discontinuity in the velocity of the 
cavity motion p2|l at the beginning of the integration does not play an important role 
for the long time behaviour of the particle production. 



Finally, we show one example for particle creation caused by the cavity motion 



to illustrate an interesting effect which appears in this case. Figure El shows the 

f Note that we still work with the same vacuum state |17o) as before (the one defined with respect to 
the eigenfrequencies fij^) and do not change the notion of particles. Instead we change only the external 
frequency Wcav which now is not twice the frequency of a quantum mode defined with respect to the 
vacuum state |0o). 



44- m = 5s{t) 




(29) 
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Figure 17. (a) Number of particles produced in a cavity vibrating with H28|l but 
cavity frequency u>cav = '^j^jy+T^j i-S- no detuning, for parameters = 1, e = 0.001 
and fcmax = 30, together with the analytical predictions H23|l and 124|l . (b) Particle 
spectrum corresponding to (a) for the two cut-off parameters femax — 20 and fcmax — 30. 




Figure 18. (a) Number of particles produced in a cavity vibrating with 1)28(1 
but cavity frequency Wcav = 2 ^^^j^^^^ , i.e. no detuning, for parameters Iq = 1, 
e = 0.01 and fcmax — HO, together with the analytical predictions H23II and H24|l . 
(b) Particle spectrum corresponding to (a) for the two cut-off parameters fcmax = 100 
and fcmax = 110. 



numerical results obtained for the parameters /q = 1, e = 0.001 and tUcav = The 
qualitative behaviour of the time evolution of the number of created particles looks 
similar to the one observed for the cavity motion (j2HI)- But in the case of the motion 
(j29|l the shape of the particle spectrum changes drastically compared to the scenarios 
which we have discussed before [see Fig. UHl (b)]. Now, not only odd but also even 
modes become excited. The resonance mode is of course the one which is excited most 
but a kind of regular structure in the spectrum seems to be absent. It is suggested that 
the reason for this irregular spectrum is due to the more complicated structure of the 
coupling function 7(t) = for the motion compared to and (f2H|) [see Fig. [T] 
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Figure 19. (a) Number of particles produced in a cavity vibrating with (|29|) and 
cavity frequency Wcav = for parameters Iq — I, t = 0.001 and femax ~ 50. (b) 
Particle spectrum corresponding to (a) for the two cut-off parameters fcmax = 40 and 



(b)]. To investigate such phenomena will be part of further work. 



5. Conclusions 



We have presented a parametrization for the time evolution of the field modes inside 
a dynamical cavity allowing for efficient numerical calculation of particle production in 
the dynamical Casimir effect. The creation of real massless scalar particles in a one- 
dimensional vibrating empty cavity has been studied numerically for three particular 
wall motions taking the intermode coupling into account. The comparison of our 
numerical results obtained for the cavity motion (j22I) with the analytical predictions 
derived in ^1 ^] shows that the numerical calculations are reliable and that 
the introduced method is appropriate to study particle creation in dynamical cavities. 
Furthermore, the range of validity of the analytical expressions found in J3j has been 
investigated. We have found that the Eqs. p3|l and (j24|l derived in the limit e -C 1 hold 
to describe the numerical results up to vibration amplitudes of the order of e = 0.01. 
The impact of the discontinuity in the velocity of the mirror at the beginning of the 
motion yielding instantaneous particle creation has been discussed. These results have 
been compared with the results obtained for the cavity motion ()28p . The numerical 
simulations show that the short time behaviours of the number of created particles 
are similar in both cases, i.e. N{t) grows quadratically with time. But the long 
time behaviours are different. While for the cavity motion ()22|) the total particle 
number increases quadratically for large times (r > 1) it passes trough a linear regime 
(i < r < 1) for the motion ()28p and afterwards slows down further. The different 
behaviours of the particle creation for the two cavity motions (j22j) and (j28p yield a 
different behaviour for the time evolution of the energy associated with the created 
quantum radiation. We have found that, in accordance with [15J, the energy in a 
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cavity vibrating with (|22|) increases exponentially whereas, in the time range under 
consideration, it obeys a quadratic law (for e = 0.001 and 0.01) in the case of the 
cavity motion ()28j) . The explanation for the differences in the behaviour of the particle 
production for the cavity motions and is that, due to detuning, Eq. (PHj) is not 
a real resonant cavity motion. Without detuning, both cavity motions yield the same 
qualitative behaviour for the time evolution of the particle number. Taking this into 
account, we may conclude that the discontinuity in the velocity of the cavity motion 
(I22|l at the beginning of the dynamics does not strongly affect the long time behaviour 
of the particle production. 

An extension of the presented method to massive fields as well as to higher 
dimensions is straightforward, which makes it possible to study the particle creation in 
the dynamical Casimir effect for scenarios where no analytical results can be deduced. 
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